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Abstract 

We report the crystal structure of a new polymorph of Ca3P2, and an analysis of its electronic 
structure. The crystal structure was determined through Rietveld refinements of powder syn¬ 
chrotron x-ray diffraction data. Ca3P2 is found to be a variant of the Mn5Si3 structure type, with 
a Ca ion deficiency compared to the ideal 5:3 stoichiometry to yield a charge-balanced compound. 
We also report the observation of a secondary phase, Ca5P3H, in which the Ca and P sites are 
fully occupied and the presence of interstitial hydride ions creates a closed-shell electron-precise 
compound. We show via electronic structure calculations of Ca3P2 that the compound is stabilized 
by a gap in the density of states compared to the hypothetical compound Ca5P3. Moreover, the 
calculated band structure of Ca3P2 indicates that it should be a three-dimensional Dirac semimetal 
with a highly unusual ring of Dirac nodes at the Fermi level. The Dirac states are protected against 
gap opening by a mirror plane in a manner analogous to graphene. The results suggest that further 
study of the electronic properties of Ca3P2 will be of interest. 
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INTRODUCTION 


Through studies of the Ca-P chemical system, the compounds CaP, Ca3P2, CaPs, CaPs, 
and CasPs have been reported to exist, but the full phase diagram is not known P^G]. In 
particular, knowledge of the compound Ca3P2 is lacking. Ca3P2 is produced in industrial 
quantities and has been used in several applications, including as a rodenticide [ 7 ], an in¬ 
gredient in pyrotechnic bombs | 8 ] , and a synthetic precursor to other phosphide compounds 
[ 0 ]. Computational studies have assumed that it crystallizes in the same structure type 
as Mg 3 P 2 ra. but to the best of the authors’ knowledge, the full crystal structure Ca3P2 
has never been reported. Furthermore, although the electronic characteristics of numerous 
semiconducting metal phosphide compounds such as InP [mill], GaP [l 3 l[Tl], and Cd3P2 
[ 15 ], have been characterized extensively, to date, there has been limited investigation of 
the electronic properties of calcium phosphides. The current study is motivated by the non¬ 
existence of a complete phase diagram for the Ca-P chemical system as well as the gap in 
our knowledge of the properties of calcium phosphides. 


Here, we report the existence, crystal structure, and electronic structure analysis of a 
polymorph of Ca3P2 that crystallizes in the Mn5Si3 structure type (space group PGs/mcm), 
distinct from a previously reported structure for Ca3P2 | 3 |. Rietveld refinements of high- 
resolution synchrotron x-ray diffraction data indicate that incomplete occupancy of the 
Ca sites yields the Ca^^-P^“ charge-balanced compound. We also report a new minority 
phase in the same space group corresponding to the formula Ca5P3H, in which interstitial 
hydride ions stabilize the full occupancy of the Ca sites by balancing the overall charge of the 
compound. Finally, we present electronic structure calculations predicting that Ca3P2 in the 
Mn5Si3 structure type is a three-dimensional ( 3 D) Dirac semimetal with a highly unusual 
ring of electronic states with Dirac (i.e. linear E vs. k) dispersions crossing the Fermi energy. 
This Dirac ring, comprising states in the k^-ky plane, is protected by crystalline reflection 
symmetry and negligibly gapped by spin-orbit coupling in a manner highly analogous to the 
Dirac states of graphene usmu. Ca3P2 may therefore be a promising candidate for the 
further study of novel quantum transport phenomena in 3 D systems. 
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EXPERIMENTAL 


Calcium metal pieces (99.5% purity) and red phosphorus powder (98.9% purity) were 
combined in a stoichiometric ratio of 5:3 and sealed in a welded tantalum tube, with no 
additional purihcation measures undertaken. The tube was heated in an induction furnace 
for 2 hours at 1200 °C as measured by a pyrometer. The resulting product was a hne blue- 
black powder that was sensitive to moisture. Commercial grade Ca 3 P 2 was obtained from 
Alfa Aesar (CAS 1305-99-3) for comparison purposes. Initial structural investigation was 
carried out via powder x-ray diffraction (XRD) on a Bruker D8 Focus diffractometer with 
Cu Ka radiation. Samples were placed in an airtight holder. 

High-resolution synchrotron XRD data was collected using the IIBM beamline at the 
Advanced Photon Source (APS) of Argonne National Laboratory using a wavelength of 
0.41384 A and a data collection range of 26 = 0.5 degrees to 50 degrees. Samples were hrst 
loaded into 0.5 mm diameter glass capillaries and sealed in an argon atmosphere due to 
their air-sensitive nature. The capillaries were placed in 0.8 mm diameter Kapton tubes for 
measurement. The obtained XRD patterns were rehned using the Rietveld method [IH] and 
the program FullProf |l9j . 

Electronic structure calculations were performed in the framework of density functional 
theory (DFT) using the wien2k [20] code with a full-potential linearized augmented plane- 
wave and local orbitals [FP-LAPW -|- lo] basis [2TH23] together with the Perdew-Becke- 
Ernzerhof (PBE) parameterization [2l| of the Generalized Approximation (GGA) as the 
exchange-correlation functional. The plane wave cut-off parameter Rmt-^max was set to 7 
and the irreducible Brillouin zone was sampled by 133 fc-points. The partial occupancy was 
modeled with the Virtual Grystal Approximation (VGA). 

Grystal orbital Hamilton population (GOHP) calculations was calculated via Tight- 
Binding Linear-Mufhn-Tin-Orbital Atomic Sphere Approximation (TB-LMTO-ASA) using 
the Stuttgart code |25l |26|. A mesh of 6 x 6 x 4 fc-points in the irreducible wedge of the hrst 
Brillouin zone was used to obtain GOHP curves. 
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RESULTS AND DISCUSSION 


Structural refinement and composition of Ca3P2 and CasPsH 

Preliminary Rietveld refinements of powder XRD patterns taken on a laboratory diffrac¬ 
tometer suggested that the obtained compound crystallized in the space group PGs/mcm, 
analogous to the compound CasAss m These diffraction patterns were distinct from that 
of commercially available Ca 3 P 2 , which was amorphous. 

Rietveld rehnements of the high-resolution synchrotron powder XRD patterns were car¬ 
ried out in the space group PQ'i/mcm (Table 1^ Table [n| and Figure [2H]. Close inspection 
of the hOO and 00/ peaks revealed anisotropy (which could not be explained by symmetry¬ 
lowering) consistent with a two-phase model (see inset of Figure [^. Note that this small shift 
in lattice parameters could not be seen on a lab diffractometer and requires high-resolution 
diffraction data. This second, minority phase occurs in the same space group as the majority 
phase with slightly different a and c lattice constants and was necessary to accurately model 
the pattern (inclusion of a second phase decreased the overall from approximately 22 to 
less than 15). Further systematic peak broadening of the 00/ reflections was modeled with 
selective peak broadening in the last stages of the rehnement. 

The freely rehned site occupancies of Ca and P for the majority and minority phases 
converged to ratios within error of 3:2 and 5:3, respectively. The majority phase contributes 
roughly 80% by weight, whereas the minority phase is approximately 20%. Based on Corbett 
and Leon-Escamilla’s previous characterization of stuffed hydride Mn 5 Si 3 -type compounds, 
such as CasAssH [29l [30], the minority phase was modeled as CasPsH. We attribute the 
presence of interstitial H to hydride impurities in the reactant alkaline earth metal or to H 2 
diffusion through the Ta tubes at high temperature, both of which are common sources of 
contamination [29|. While we acknowledge that x-ray diffraction is expected to be insensitive 
to H-scattering, the inclusion of H on the 2b site is necessary through chemical reasoning 
to obtain a charge-balanced compound for the minority phase, discussed in further detail 
below. 

The minority phase was rehned with an overall thermal parameter, whereas the atoms 
in the majority phase were treated with isotropic thermal parameters (the Ca sites were 
constrained to have the same thermal parameter). In the majority phase, the composition 
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was determined by performing a series of refinements where the occupancies of the Ca 
sites (constrained to have the same percent occupancy) were hxed at incrementally changed 
values, and the other parameters were allowed to freely rehne. This accommodates the 
possible presence of correlations between rehnement parameters EH. As shown in the inset 
of Figure the -Rsragg for the majority phase reached a minimum when the Ca sites were 
90-95% occupied, yielding a formula within error of Ca 3 P 2 . The occupancies of the Ca 
sites were therefore hxed to stoichiometric Ca 3 P 2 in the hnal stages of the rehnement. No 
evidence of vacancy ordering or P dehciency was found. The minority phase consistently 
rehned to full occupancy of the Ca and P sites (the occupancy of H was not rehned), yielding 
the charge-balanced compound Ca 5 P 3 H. Additional tests were performed with oxygen as the 
interstitial ion stabilizing the minority phase, but this led to poorer agreement of the model 
to the data, and thus that possibility was eliminated. The hnal atomic coordinates and 
occupancies for both phases are given in Table Two impurity phases are also present in 
small proportions, Ca4P20 (2.5%) and Ca^jTaOs (0.02%). No attempt was made to obtain 
the minority hydride phase in pure form. 

Visualizations of the structures of Ca 3 P 2 and CasPsH are shown in Figure]^ The lattice 
parameters of the phases are similar (Table |^, with Ca 5 P 3 H having a slightly smaller unit 
cell volume (0.1%). This is in line with the volume trends observed by Leon-Escamilla 
and Corbett for alkaline earth arsenides of the MnsSis-type structure in the presence of 
interstitial hydrogen 1301 , albeit the volume change here is to a lesser degree. We attribute 
this to the fact that the Ca dehciency in the hydrogen-free compound reduces the volume 
of this phase relative to the other alkaline earth pnictides with all sites fully occupied. 

As Leon-Escamilla and Corbett have demonstrated, the inclusion of hydrogen in Mn 5 Si 3 - 
type phases serves a chemically straightforward stabilizing function [291130] . Electron count¬ 
ing reveals that a compound with the formula A 5 B 3 , where A is an alkaline earth metal and 
5 is a pnictogen, has one extra valence electron per formula unit (assuming the oxidation 
states of +2 for A and -3 for B). The inclusion of an interstitial atom such as H or F 
adopting a -1 oxidation state yields a closed-shell, charge-balanced system. For the light 
and relatively polarizing elements Ca and P, we expect a metallic CasPs phase to be more 
unstable than the analogous A^B^ phases with heavier elements. The above electron count¬ 
ing scheme can be applied to qualitatively rationalize both Ca dehciency in the Ca 3 P 2 phase 
and the H interstitial in the minority Ca 5 P 3 H phase. 
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Previous thermodynamic studies of the Ca-P chemical system have noted that systematic 
direct study of the phase behavior is difficult owing to the high volatility and reactivity of 
the constituents |2]. Based on conditions for Ca 3 P 2 synthesis given in the existing literature 
PE], we hypothesize that the MnsSis-type compound reported here may be a higher- 
temperature phase favored by the sealed environment of the tantalum reaction container. 
Detailed structural characterization of the polymorph reported in previous studies |3] as well 
as additional examination of the thermodynamic properties will be necessary to explore the 
exact relationship between these phases and others in the Ca-P system. 


Electronic structure of Ca3P2 

To further investigate the nature of stability and bonding in Ca 3 P 2 , quantitative calcula¬ 
tions of the electronic band structure, density of states (DOS), and crystal orbital Hamilton 
population (COHP) were carried out (Figure [^. Both hypothetical Ca 5 P 3 and Ca 3 P 2 were 
calculated in order to confirm the validity of the electronic structure model for this system. 

Both the COHP and the DOS indicate a region of electron count stability for Ca 3 P 2 : the 
Fermi level lies in a deep pseudogap. The Fermi level of a hypothetical Ca 5 P 3 compound 
would be located near a peak in the DOS, indicating chemical instability (Figure]^). Anal¬ 
ysis of the COHP indicates that only bonding states are occupied in Ca 3 P 2 (Figure]^). 
The additional filled states in CasPs are bonding in nature as well, but the larger COHP 
value indicates the composition is less stable. 

Figure and Figure show the band structure of hypothetical Ca 5 P 3 and Ca 3 P 2 
calculated with the VCA. The features of the electronic structures are shifted in energy but 
otherwise analogous, showing that a rigid band model is valid in this case. In Ca 3 P 2 , two 
linear crossings of bands are observed at the Fermi energy along the the F-K and the 
F-M lines. In hypothetical Ca 5 P 3 , the same feature is present 0.5 eV below E-p. 

In the PQ^/mcm space group, crystalline symmetry protects these crossings. Along the 
F-K and F-M lines, the valence band has irreducible representation Fi and the conduction 
band has representation F 2 . Due to these different irreducible representations, the band 
crossings are protected by a mirror plane and C 2 rotational symmetry in the absence of 
spin-orbit coupling (SOC). SOC allows mixing of the states and opens the crossing [32]. 
However, the effect of SOC is negligibly small for the light elements Ca and P. Therefore, 
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even when SOC is applied to our calculation, we do not observe a gap in the crossing. 

This type of linear band crossing, also called a bulk Dirac cone, at the Fermi energy 
has been of recent interest due to interesting physical properties observed in compounds ex¬ 
hibiting such a band structure, such as high magnetoresistance and extremely high mobility 
[33H38] . Graphene, for example, is a 2D Dirac semimetal [HI ESI HQ]. One of the prominent 
examples of a 3D Dirac semimetal is Cd 3 As 2 , which has recently gained a lot of attention 
ESHSZl EHH!. Although Cd 3 As 2 crystallizes in a very different crystal structure [H], the 
chemical relation to Ca 3 P 2 is still visible. Ca 3 P 2 may be a less toxic alternative to Cd 3 As 2 . 
The other known 3D Dirac semimetal, Na 3 Bi, has the disadvantage of being extremely air 
sensitive, which makes handling very difficult [38]; Ca 3 P 2 is air-sensitive as well, but to a 
much lesser extent. 

Figure shows a three-dimensional plot of the band structure of Ca 3 P 2 in the and 
ky directions vs. E. The band crossings with linear dispersions along F-K and F-M form 
a closed ring around the Cq rotation axis at the Fermi energy. As a result of the C 21 ; 
symmetry along these lines and the different irreducible representations of the bands, the 
Dirac crossings are symmetry-protected (without SOC) by the k^ = 0 mirror plane and 
C 2 rotational symmetry. For an arbitrary line connecting F-K and F-M, the point group 
is Cg and the mirror plane remains a protecting element for all Dirac nodes within the 
kx-ky plane. Ca 3 P 2 is therefore the hrst 3D Dirac semimetal predicted where a mirror 
plane is the symmetry element that protects the crossing of the bands. Mirror planes have 
been suggested as a protecting symmetry element by Chiu and Schnyder, who propose a 
classihcation scheme for reflection-symmetric Dirac semimetals, among which graphene is an 
example [TB]. An analogy therefore can be drawn between the band structures of graphene 
and Ca 3 P 2 : In graphene, six discrete Dirac points are present at the K and K’ points within 
the kx-ky plane [171 EHl SO]. In Ca 3 P 2 , a continuous, two-dimensional ring of Dirac nodes 
with sixfold rotational symmetry exists in the kx-ky plane. Figure and c illustrate the 
relationship. 

Additionally, we calculated the band structure of several other compounds in the Mn 5 Si 3 
structure type and found the same kind of Dirac crossings 0.5 eV below the Fermi level, 
similar to hypothetical Ca 5 P 3 . The vacancies in Ca 3 P 2 in the Mn 5 Si 3 structure allow for 
the Dirac states to be at the Fermi level. A hypothetical compound such as Ca 4 NaP 3 , for 
example, should also have Dirac states at the Fermi level assuming a rigid band model. 
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Modulating the electron count of other compounds in this structure type, especially those 
with light elements, may be a pathway to new 3 D Dirac semimetals. Other reported topo¬ 
logical materials such as Ru2Sn3 [ 3 S] and [Tl4](Tli_3,Sn2:)Te3 [IB] also have in-plane mirror 
planes causing a Dirac crossing if SOC is not applied in the calculation. The large SOC in 
these materials gaps out these crossings, however. These results imply that in-plane mirror 
planes can generally produce new 3 D Dirac materials if SOC is small enough. 


CONCLUSIONS 

The results of the structural rehnements suggest that Ca vacancies in a Ca-P compound 
in the Mn5Si3 structure type yields a composition of Ca3P2, with the driving force being the 
greater stability of a charge-balanced, closed-shell composition relative to the hypothetical 
compound Ca5P3. Electronic structure calculations further support the stability of Ca3P2 
compared to Ca5P3. The hndings suggest that further examination of the equilibrium com¬ 
position of compounds in the Mn5Si3 structure type, such as other alkaline earth metal 
phosphides and arsenides, may be of future interest. In addition, study of the exact mech¬ 
anism and extent of hydrogen interstitials in stuffed Mn5Si3-type structures, drawing upon 
the existing work of Leon-Escamilla and Corbett | 29 lEQ], may provide additional insight 
into the composition and stability of these phases and their parent compounds. 

The calculated band structure of Ca3P2 suggests the compound may be a candidate for 
realizing a 3 D Dirac semimetal ground state. To the authors’ knowledge, this is the hrst 
proposed Dirac semimetal in this structure type, as well as the hrst 3 D compound predicted 
to have bulk Dirac states protected by crystalline rehection symmetry that are negligibly 
gapped by SOC. The similarities between Ca3P2 and graphene, in particular, merit further 
study of its electronic properties. The result illustrates that rehection symmetry can operate 
as a protecting element for realizing 3 D Dirac semimetals, which we hope will motivate 
additional exploration of such systems. 
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FIG. 1. Rietveld refinement of synchrotron XRD data. Peaks corresponding to the phases Ca 3 P 2 
and CasPsH, as well as Ca4P20 and Caa^TaOa impurity phases are indicated. Left inset: ^Bragg 
values for the majority phase refinement vs. Ca occupancy, wherein the lowest values of ^Bragg 
correspond to Ca occupancies between 90-95%. Right inset: examples of peak splitting fitted by a 
two-phase model. 
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TABLE I. Data of crystallographic refinement for Ca 3 P 2 and CasPaH. f Sum of weight fractions 
is slightly less than 100 % due to the presence of small amounts of impurity phases. 


Crystal system 

Hexagonal 

Overall refinement 

Space group 

P 63 / mcm 


Temperature (K) 

295.0 


Z 

2 


Rwp 

25.0 


Rp 

23.1 



12.8 



Ca3P2 

CasPsH 

Rf 

11.1 

11.9 

R'Bragg 

12.3 

13.8 

a (A) 

8.25593(1) 

8.25979(2) 

c (A) 

6.83551(2) 

6.82308(3) 

V (A3) 

403.491(1) 

403.134(2) 

Density (g/cm^) 

2.249 

2.425 

Weight fractionf (%) 

79.1(8) 

18.4(2) 
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TABLE II. Atomic coordinates and parameters in the crystal structure of Ca 3 P 2 and CasPsH. 


Atom 

Site 

X 

y 

z 

B 

Occ. 

Phase: Ca 3 P 2 





Biso 


Cal 

4d 

1/3 

2/3 

0 

1.62(3) 

1 

Ca2 

65 

0.2544(2) 

0 

1/4 

1.62(3) 

0.9 

P 

65 

0.6096(2) 

0 

1/4 

0.56(3) 

0.9 

Phase: CasPsH 




-Bov 


Cal 

4d 

1/3 

2/3 

0 


1 

Ca2 

6g 

0.2543(3) 

0 

1/4 

2.2(4) 

1 

P 

6g 

0.6104(1) 

0 

1/4 


1 

H 

2b 

0 

0 

0 


1 
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FIG. 2. The crystal structures of Ca 3 P 2 and CasPsH. a. View of Ca 3 P 2 along the [001] direction, 
b. Unit cell of Ca 3 P 2 with planar network of Ca-P atoms indicated, c. Representation of H 
interstitial atoms in Ca 5 P 3 H in octahedral coordination with Ca atoms. 
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FIG. 3. a. and b. The density of states (DOS) and the crystal orbital Hamilton population 
(COHP) of Ca 3 P 2 calculated with the VGA. The Fermi level of hypothetical CasPs (with all Ga 
sites fully occupied) is indicated, c. and d. The calculated electronic band structures of CasPs (c.) 
and Ca 3 P 2 (d.), showing bulk Dirac cones between F and M, and F and K. Each color corresponds 
to a different irreducible representation, e. A diagram of the band inversion at the Fermi energy 
in Ca 3 P 2 . 
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FIG. 4. a. Three-dimensional plot of the band structure of Ca 3 P 2 along kx and ky vs. E. The 
closed ring formed by the band crossings in the kx-ky plane at the Fermi energy is highlighted in 
red. b. and c. Graphical representation of the band crossings in the kx-ky plane of the Brillouin 
zone for Ga 3 P 2 (b.) and graphene (c.). 
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